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ABSTRACT.  We  study  plane  strain  thermomechanical  deformations  of  a  prismatic 
viscoplastic  body  of  square  cross-section  and  deformed  at  a  nominal  strain-rate  of  5000 
sec'1 .  The  body  has  two  thin  layers  placed  symmetrically  about  the  horizontal  centroidal 
axis.  The  layer  material  differs  from  that  of  the  body  in  only  the  value  of  the  yield  stress  in 
a  quasistatic  simple  compression  test.  The  yield  stress  for  the  layer  material  is  taken  to  be 
either  one-fifth  or  five  times  that  of  the  matrix  material.  Three  cases,  namely,  when  there  is 
an  elliptical  void  with  <ts  major  axis  aligned  along  the  horizontal  centroidal  axis  or  the  ver¬ 
tical  centroidal  axis,  but  the  void  center  coincides  with  the  center  of  the  cross-section,  and 
when  there  are  two  elliptical  voids  with  major  axes  aligned  along  the  vertical  centroidal  axis 
and  a  void  tip  abuts  the  layer/matrix  interface  are  studied.  The  deformations  are  assumed 
to  be  symmetrical  about  the  vertical  and  horizontal  centroidal  axes. 

It  is  found  that  in  each  case  shear  bands  initiate  from  points  on  the  traction  free  edges 
where  the  matrix/layer  interfaces  intersect  them  and  propagate  into  the  softer  material.  For 
the  soft  layer  these  bands  initially  merge  into  one  and  propagate  horizontally. 

Subsequently,  each  of  these  bands  bifurcates  into  two  which  propagate  into  the  matrix 
material  along  the  direction  of  the  maximum  shear  stress.  There  is  minimal  interaction 
between  these  bands  and  those  initiating  from  points  near  the  void  tips. 

INTRODUCTION.  Adiabatic  shear  bands  are  narrow  regions  of  intense  plastic 
deformation  that  form  during  high  strain-rate  processes,  such  as  shock  loading,  ballistic 
penetration,  metal  forming,  and  machining.  As  these  bands  generally  precede  material 
fracture,  a  knowledge  of  factors  that  inhibit  or  enhance  their  giowtli  is  essential  to  the  pro¬ 
duction  of  durable  materials  and  more  efficient  manufacturing  processes.  These  bands 
form  in  both  ferrous  and  nonferrous  alloys. 

Johnson  (1987)  has  recently  pointed  out  that  Tresca  (1878)  and  Massey  (1921) 
observed  hot  lines,  now  referred  to  as  adiabatic  shear  bands,  during  the  forging  of  platinum. 
Both  Tresca  and  Massey  stated  that  these  were  the  lines  of  greatest  sliding,  and  also  there¬ 
fore  the  zones  of  greatest  development  of  heat.  Wulf  (1978)  has  reported  experimental 
observations  of  adiabatic  shear  bands  in  high  strain  rate  (2000  to  25000  sec'1)  compression 
of  7039  aluminum  armour.  He  found  that  the  cross-section  of  the  cylindrical  specimens 
changed  from  circular  to  elliptical  after  the  compression  test,  and  adiabatic  shear  bands 
formed  in  the  specimens  which  subsequently  failed  by  crack  propagation  along  the  dom¬ 
inant  band.  Further  references  to  the  analytical,  numerical  and  experimental  work  on  shear 
banding  may  be  found  in  two  recent  books  (Dodd  and  Bai  (1987),  Semiatin  and  Jonas 
(1984)). 

Recently,  LeMonds  and  Needleman  (1986),  Needleman  (1989),  Anand  et  al.  (1988), 
Zbib  and  Aifantis  (1988),  Batra  and  Liu  (1989,1990),  Zhu  and  Batra  (1990),  Batra  and  Zhu 
(1990),  and  Batra  and  Zhang  (1990)  have  studied  the  phenomenon  of  shear  banding  in 
)  plane  strain  deformations  of  a  viscoplastic  solid.  Whereas  Needleman  studied  a  purely 

mechanical  problem,  other  works  have  treated  a  coupled  thermomechanical  problem. 
LeMonds  and  Needleman,  Zbib  and  Aifantis,  and  Anand  et  al.  neglected  the  effect  of  iner¬ 
tia  forces  on  the  ensuing  deformations  of  the  body.  These  investigations  have  employed 

Supn  ir  teo  oy  die  u.  S.  Army  Research  Office  Contract  DAAL03-88-K-0184  to  the  Uni¬ 
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different  constitutive  relations,  different  technique-  to  integrate  the  stiff  set  of  governing 
partial  differential  equations,  and  have  generally  assumed  that  the  entire  body  or  the  por¬ 
tion  of  the  body  whose  deformations  were  analyzed  had  only  one  defect  in  it.  The  prismatic 
body  whose  plane  strain  thermomechanical  deformations  are  studied  herein  is  of  a  square 
cross-section  and  has  two  thin  layers  made  of  a  viscoplastic  material  different  from  that  of 
the  body  and  placed  symmetrically  about  and  parallel  to  the  horizontal  centroidal  axis. 
These  horizontal  layers  may  be  thought  of  as  representing  planes  of  chemical  inhomoge¬ 
neity.  Also,  as  stated  above  in  the  abstract,  there  is  either  an  elliptical  void  at  the  center  of 
the  cross-section  or  two  ellipsoidal  voids  with  major  axes  along  the  vertical  centroidal  axis 
and  tips  touching  the  layer/matrix  interfaces.  The  voids  can  form  during  manufacturing 
However,  the  symmetrical  situation  considered  herein  is  to  simplify  the  problem.  The  ver¬ 
tices  of  the  ellipsoidal  void  on  its  major  axis  and  the  points  on  the  free  edges  where  the  thin 
layer  and  the  matrix  materials  meet  act  as  nuclei  for  the  initiation  of  shear  bands.  It  thus 
becomes  an  interesting  exercise  to  study  the  initiation  and  growth  of  various  bands  and  the 
interaction,  if  any,  amongst  them.  We  account  for  the  effect  of  inertia  forces,  strain-rate 
sensitivity  of  the  materials,  thermal  softening  effects,  heat  conduction,  and  the  heat  gener¬ 
ated  due  to  plastic  working. 


FORMULATION  OF  THE  PROBLEM.  Figure  1  depicts  the  cross-section  of  the 
prismatic  body,  and  the  relative  dimensions  of  the  ellipsoidal  void  and  the  two  thin  layers. 
For  this  case,  the  centers  of  the  void  and  the  cross-section  coincide  and  the  major  axis  of  the 
void  coincides  with  the  vertical  centroidal  axis  of  the  cross-section.  It  is  assumed  that  the 
body  is  loaded  along  the  vertical  axis,  plane  strain  state  of  deformation  prevails,  and  that  the 
deformations  are  symmetrical  about  the  two  centroidal  axes.  Thus  the  deformations  of  the 
material  in  the  first  quadrant  are  analyzed.  We  use  a  fixed  set  of  rectangular  cartesian  coor¬ 
dinate  axes  and  the  referential  description  of  motion  to  describe  the  thermomechanical 
deformations  of  the  body.  The  governing  equations  are: 


balance  of  mass:  (pJ) '  =  0,  (2.1) 

balance  of  linear  momentum:  p0vj  =  Tjaa,  (2.2) 

balance  of  moment  of  momentum:  Tjtaxj  a  =  Tj>axja,  (2.3) 

balance  of  internal  energy:  p0e  =  -Qace  +  Tjavja,  (2.4) 

where 

Tia  =  (Pq/p)  aij  xa,j/ 

°  ij  =  -  B  (p/p0  —  lj  S  2  j  +  2p  j  i  (2.5) 

2U  =  [a0/j3I]  (1  +  bl)m  (1  -  ced),  (2.6) 

I2  =  (1/2)  Dij  Dij,  (2.7) 

Dij  =  ~  Dkk  ^  ii  >  (2.8) 

Qqs  =  (Po/P)  ,  i •  <Ji  =  ~  k  0  ,  ii  (2.9) 

e  =  c  9  +  B(p/p0  -  1)  p/p2  .  (2.10) 


In  these  equations  xj  gives  the  position  at  time  t  of  the  material  particle  X^,  vj  =  tj  is 
its  velocity  in  the  xj-direction,  p  is  its  present  mass  density,  p0  its  mass  density  in  tne  refer¬ 
ence  configuration,  J  =  det  [xj  a],  xj((J  =  d\\/d  X^,  Tja  is  the  first  Piola-Kirchoff  stress  ten¬ 
sor,  ojj  is  tne  Cauchy  stress  tensor,  e’is  the  specific  internal  energy,  Qa  is  the  heat  flux 
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measured  per  unit  area  in  the  reference  configuration,  D  is  the  strain-rate  tensor  and  tJ  is 

its  deviatoric  part,  a  superimposed  dot  indicates  material  time  derivative,  a  comma  followed 
by  index  a  (j)  implies  partial  differentiation  with  respect  to  Xa  (x;),  and  a  repeated  index 
implies  summation  over  the  range  (1,2)  of  the  index.  In  the  constitutive  relations  (2.5), 

(2.9)  and  (2.10),  the  material  parameter  B  represents  the  bulk  modulus,  a  0  is  the  yield 
stress  in  a  quasistatic  simple  compression  test,  parameters  b  and  m  characterize  the  strain- 
rate  sensitivity  of  the  material,  a  describes  its  thermal  softening,  6  equals  the  temperature 
change  of  a  material  particle  from  that  in  the  reference  configuration,  k  is  the  constant  ther¬ 
mal  conductivity  and  c  is  the  constant  specific  heat.  Here  we  have  not  considered  the 
stresses  caused  by  the  thermal  expansion. 

The  foregoing  equations  hold  in  regions  occupied  by  the  matrix  and  the  layers.  The 
values  of  material  parameters  for  the  matrix  and  the  layer  materials  are  the  same  except 
that  either 

a Q  layer  =  5  oQ  matrix/ 
or 

0O  layer  =  (1/5)  aQ  matrix. 

In  terms  of  the  deviatoric  stress  s  defined  by 

s  =  o  +  [B(p/p0  -  1)  -  (2p/3)  tr  D]l, 

2p  D , 

equations  (2.12),  (2.5)  and  (2.6)  give 

(1/2  tr  s2)V2  =  (<70/j  3)  (1  -  ad)  (1  +  bl)m. 

We  assume  that  the  body  is  initially  at  rest  at  a  uniform  temperature,  has  a  constant 
mass  density  and  is  initially  stress  free.  That  is 

P  ( x » o)  =  1,  v ( x , o)  =  o,  0(x,o)  =  0  .  (2.14) 

For  the  material  in  the  first  quadrant,  we  impose  the  following  boundary  conditions. 


v2 

=  -  h(t)  ,  T12  =  d  and 
surface  AB, 

Q2  =  0, 

on 

the 

top 

(2.15) 

Til 

=  0,  T21  =  0  and  Q1  = 
surface  BC, 

=  o, 

on 

the 

right 

(2.16) 

v2 

=  Of  T]_2  =  0  and  Q 2  = 
surface  CD, 

0, 

on 

the 

bottom 

(2.17) 

Tic* 

Nqj  =  0 ,  Qa  Na  —  0 , 

surface  DE, 

on 

the 

void 

(2.18) 

V1 

=  0,  T21  =  0  and  = 

surface  EA. 

0, 

on 

the 

left 

(2.19) 

That  is  the  top  surface  is  moving  downward  with  a  speed  h(t),  contact  between  it  and  the 
loading  device  is  smooth,  the  right  surface  is  traction  free,  and  the  entire  boundary  is  ther- 


(2 . lla) 
(2.11b) 

(2.12a) 

(2.12b) 

(2.13) 
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mally  insulated.  If  at  any  time  during  the  deformations  of  the  body,  a  point  on  the  void 
surface  touches  the  vertical  axis,  the  boundary  condition  on  it  is  changed  to  (2.19).  The 
boundary  conditions  (2.17)  and  (2.18)  reflect  the  presumed  symmetry  of  deformations 
about  the  x\  and  X2  axes.  For  the  loading  function  h(t)  we  take 
h(t)  =  v0  t/tr,  o  <;  t  i  tr, 

=  vQ  t  >  tr.  (2.20) 

The  steady  speed  v0  of  the  top  surface  of  the  block  is  reached  in  time  u. 

The  matrix  and  the  layer  are  assumed  to  be  perfectly  bonded.  Thus  at  the  common 
interfaces  between  them,  the  velocity  field,  surface  tractions,  the  temperature  and  the  nor¬ 
mal  component  of  the  heat  flux  are  assumed  to  be  continuous. 

For  other  configurations  of  the  voids,  the  boundary  conditions  are  appropriately 
modified. 

RESULTS.  The  finite  element  code  developed  by  Batra  and  Liu  (1989, 1990)  was 
modified  to  analyze  the  present  problem.  In  order  to  compute  results,  we  used  the  follow¬ 
ing  values  of  vanous  material  and  geometric  parameters. 

b  =  10000  sec,  oQ  =  333  MPa,  k  =  49.22  W  nT1  °C_1, 
m  =  0.025,  c  =  473  J  Kg"1  “C"1,  pQ  =  7860  Kg  in"3, 

B  =  128  GPa,  H  =  5  mm,  vQ  =  25  m  sec"1,  (3.1) 

a  =  0.0025  "C"1. 


Thus  the  average  applied  strain-rate  -yaVg  equals  5000  sec'1,  the  reference  temperature 

dQ  m  OqKpq c)  =  89.6°  C,  and  v  -  PpVg/op  =  0.015.  The  nondimensional  number  v  signi¬ 
fies  the  effect  of  inertia  forces  relative  to  tne  flow  stress  of  the  material.  For  the  simple 
shearing  problem,  Batra  (1988)  noted  that  the  inertia  forces  play  a  noticeable  role  when  v 
=  0.004.  Thus  for  the  present  problem,  the  inertia  forces  will  very  likely  play  a  significant 
role. 

LAYER  MATERIAL  SOFTER  THAN  THE  MATRIX  MATERIAL.  Figure  2 
depicts  the  contours  of  the  maximum  principal  logarithmic  strain  e  at  an  average  strain  of 
0.079  when  the  center  of  the  ellipsoidal  void  coincides  with  the  center  of  the  cross-section 
and  the  major  axis  of  the  void  is  aligned  along  the  vertical  centroidal  axis.  These  contours 
and  other  results  reported  by  Zhu  and  Batra  (1990b)  reveal  that  shear  bands  initiate  from 

goints  on  the  right  traction  free  edge  where  the  matrix  /layer  interfaces  intersect  it. 

because  the  layer  material  is  softer  and  its  thickness  quite  small,  these  bands  merge  into  one 
band  that  initially  propagates  horizontally  into  the  layer.  When  the  matrix  material  has 
softened  somewhat  due  to  the  rise  in  its  temperature,  the  horizontally  propagating  band 
bifurcates  into  two  bands  that  propagate  into  the  matrix  along  ±45°  directions.  The  band 
propagating  into  the  matrix  material  above  the  layer  has  more  severe  deformations  associ¬ 
ated  with  it  than  the  one  propagating  into  the  matrix  material  below  the  layer. 

The  band  initiating  from  a  point  near  the  void  tip  on  the  vertical  centroidal  axis  pro¬ 
pagates  along  a  line  that  makes  an  angle  of  45°  with  the  horizontal.  This  band  seems  to 
pass  through  the  soft  layer  rather  easily. 

Figure  3  shows  the  contours  of  the  maximum  principal  logarithmic  strain  e  at  an  aver¬ 
age  strain  of  0.0333  when  the  ellipsoidal  void  is  at  the  center  of  the  cross-section  and  the 
major  axis  of  the  void  is  aligned  along  the  horizontal  centroidal  axis.  These  plots  look  quite 
similar  to  that  for  the  case  when  the  major  axis  of  the  void  coincides  with  the  vertical  cen¬ 
troidal  axis.  For  ?  f’>rther  discussion  and  details  of  results  in  this  case,  see  Batra  and  Zhu 
(1991). 

When  the  void  tip  touches  the  matrix/layer  interface  and  the  major  axis  of  the  void 
coincides  with  the  vertical  centroidal  axis,  the  contours  of  t  plotted  in  Figure  4  at  an  aver¬ 
age  strain  rate  of  0.0175  look  quite  different  from  the  previous  two  cases.  Results  given  by 
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Batra  and  Zhu  (1990)  and  these  contours  of  e  suggest  that  a  shear  band  initiates  within  the 
matrix  surrounding  the  void  tip  near  the  matrix/layer  interface  and  propagates  into  the 
matrix  material  below  the  common  interface,  the  direction  of  propagation  being  nearly  45° 
to  the  vertical  axis.  The  shear  bands  initiating  at  points  of  intersection  of  the  matrix/layer 
interfaces  with  the  right  traction  free  surface  propagate  into  the  soft  layer  and  then  bifur¬ 
cate  into  the  matrix  material  along  lines  making  an  angle  of  approximately  45°  with  the 
vertical.  The  band  in  the  layer  near  the  upper  matrix/layer  interface  bifurcates  into  the 
matrix  prior  to  that  near  the  lower  interface.  Also  the  band  in  the  layer  near  the  upper 
matrix/layer  interface  continues  to  propagate  horizontally  into  the  layer  too,  while  that  near 
the  lower  surface  does  not. 

LAYER  MATERIAL  STRONGER  THAN  THE  MATRIX  MATERIAL.  We  first 
study  the  case  when  the  center  of  the  ellipsoidal  void  coincides  with  that  of  the  cross- 
section.  In  Figure  5,  we  have  plotted  contours  of  €  at  an  average  strain  of  0.122  when  the 
major  axis  of  the  void  is  vertical.  Now  the  bands  initiating  from  points  on  the  right  traction 
free  edge  where  the  matrix/layer  interfaces  intersect  it  propagate  into  the  matrix  material 
along  lines  making  an  angle  of  approximately  ±45®  with  the  horizontal.  Recall  that  the 
matrix  material  is  softer  than  the  layer  material.  The  quarter  of  the  square  cross-section 
studied  is  divided  into  five  subregions,  each  of  which  is  deforming  essentially  rigidly,  and 
there  is  a  shear  band  at  the  four  common  boundaries  (e.g.,  see  Zhu  and  Batra  (1990b)). 
However,  when  the  major  axis  of  the  void  is  horizontal,  contours  of  e  depicted  in  Figure  6 
at  an  average  strain  of  0.031  suggest  a  picture  different  from  the  one  when  the  major  axis  of 
the  void  was  vertical.  We  should  add  that  the  average  strains  in  the  two  cases  are  quite 
different.  Hence,  a  direct  comparison  is  not  very  meaningful. 

Figure  7  shows  contours  of  e  at  an  average  strain  of  0.057  when  a  void  tip  is  at  the 
matrix/layer  interface.  These  contours  and  other  results  given  by  Batra  and  Zhu  (1990) 
reveal  that  a  shear  band  initiating  from  the  void  tip  abutting  the  matrix/layer  interface  pro¬ 
pagates  initially  along  the  interface  and  then  into  the  matrix  material  along  a  line  making  an 
angle  of  nearly  45°  with  the  vertical.  The  shear  band  initiating  from  the  lower  void  tip  also 
propagates  into  the  matrix  material  along  a  line  making  an  angle  of  approximately  45°  with 
the  vertical.  Two  shear  bands  also  initiate  from  points  on  the  right  traction  free  edge  where 
matrix/layer  interfaces  intersect  it,  and  these  bands  propagate  into  the  matrix  material  along 
lines  making  an  angle  of  45°  with  the  vertical.  Even  though  it  seems  that  near  the  vertical 
centroidal  axis  a  shear  band  has  propagated  into  the  layer,  there  is  no  localization  of  defor¬ 
mation  occurring  in  the  layer  material.  This  is  evidenced  by  the  plots  of  e  versus  the  aver¬ 
age  strain  at  several  points  in  the  layer  that  are  given  in  Fig.  8c  of  Batra  and  Zhu’s  (1990) 
paper.  Tne  contours  of  the  temperature  rise,  not  included  herein,  support  the  picture  laid 
out  above  for  the  development  of  four  bands,  two  from  void  tips  and  two  from  points  on  the 
right  traction  free  surface  where  the  layer  and  the  matrix  materials  meet. 

CONCLUSIONS.  We  have  analyzed  the  problem  of  the  initiation  and  growth  of  shear 
bands  in  a  prismatic  viscoplastic  body  containing  an  ellipsoidal  void,  and  two  thin  layers 
made  of  a  different  viscoplastic  material  placed  symmetrically  about  the  horizontal  cen¬ 
troidal  axis.  The  body  is  deformed  in  plane  strain  compression  along  the  vertical  axis  at  an 
average  strain-rate  of  5000  sec'1 ,  and  its  deformations  are  assumed  to  be  symmetrical  about 
the  two  centroidal  axes. 

Two  shear  bands  initiate  from  points  on  the  right  traction  free  edge  where  the  mat¬ 
rix/layer  interfaces  meet  it.  These  bands  propagate  into  the  softer  material.  When  the 
matrix  material  is  softer,  these  bands  propagate  along  lines  that  make  an  angle  of 
approximately  ±45°  with  the  horizontal.  However,  when  the  layer  material  is  softer,  these  bands 
essentially  merge  into  one  and  initially  propagate  horizontally  into  the  layer.  Subsequently, 
this  band  bifurcates  into  two  bands  that  propagate  into  the  matrix  material  along  ±45° 
directions. 

Shear  bands  also  initiate  from  points  near  the  void  tips  and  propagate  in  the  ±45° 
directions.  When  the  layer  material  is  stronger  than  the  matrix  material,  these  bands  do  not 
pass  through  the  layer  and  are  deflected  back  upon  arriving  at  the  matrix/layer  interface. 
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However,  they  pass  through  easily  through  a  softer  layer. 
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Fig.  2.  Contours  of  Uie  maximum  principal  logarithmic  strain  at  7avg  =  0.079.  o  „  Iayer  = 

( 1/5)  o  0  The  major  axis  of  the  ellipsoidal  void  cc  incites  with  the  verti¬ 

cal  centroidal  axis. 
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Fig.  3.  Contours  of  the  maximum  principal  logarithmic  strain  at  7ave  =  0.0333. 

0  o  layer  =  ( 1/5)  a  o  matrix,  the  major  axis  of  the  ellipsoidal  void  coincides  with 
the  horizontal  centroidal  axis. 
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